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\^ . Abstract 

m 

00 I In this paper we derive the Laplace transforms of the integral func- 

O . tionals 

'NT 

o 



p ( eMBD + 1 ) ' + q( exp(i?r) + 1 ) ] dt, 



and 



. , 1 P ( exp(i2f - 1 ) +q ( exp(i?f ) - 1 ) ) dt, 

H , 

■ where p and q are real numbers, {B^^^ : t > 0} is a Brownian motion 

with drift ^ > 0, BM(^), and {Rl ' : t > 0} is a 3-dimensional Bessel 
process, BES(3). The transforms are given in terms of Gauss' hyper- 
geometric functions and it is seen that the results are closely related 
to some functionals of Jacobi diff'usions. This work generalizes and 
completes some results of Donati-Martin and Yor jlj and Salminen 
and Yor [HI. 
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1 Introduction 

Let B^f^^ = {Bf^^ := Bt + fit : t > 0} be a Brownian motion with drift 
fi > started from 0. For the non-negative and locally integrable function / 
consider the functional 

/■oo 

Joo(/):= / fiBi'^^)ds. 
Jo 

It is known, see Engelbert and Senf jH) and Salminen and Yor [T8] . 

/oo 
f{x) dx < oo. (1.1) 

In [in] also other properties of /oo(/) are studied; in particular, existence of 
moments. 

During the recent years and because of applications in financial mathe- 
matics much interest has been focused on exponential functionals of B^f^\ We 
refer to Yor P^Hj for a collection of papers on this topic and also for further 
references. A particular role in these investigations has been played by the 
functional 

/■oo 

/ exp (-25^'^)) ds (1.2) 
Jo 

which, as shown in Dufresne 0, can be obtained as a perpetuity in a discrete 
pension funding scheme after a limiting procedure. See also Yor ^T7\ . 
In [TT] and [12] the functional 

roo 

/ {a + exp{Bi''^))-^ds (1.3) 
Jo 

is analyzed the motivation being that the functional in (jl.3ll has all the 
moments (in fact, some exponential moments) when Dufresne's functional 
()1.2p has only some moments. The distribution of the functional in ()1.3|) is 
computed in |11| for /i = 1/2 on one hand by connecting the functional to a 
hitting time of another Brownian motion with drift via random time change 
and on the other hand using the Feynman-Kac formula. 
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In this paper we derive the Laplace transform of the functional in 
for arbitrary yU > 0. In fact, we study the more general functional 




p(exp(5i^^) + 1) Vg(exp(fi;^^) + 1) dt. (1.4) 



Notice that choosing here, e.g. , p = — g = 1 we obtain the functional 



The key observation is that the functional in ()1.4|1 can be expressed, roughly 
speaking, via the first hitting time of a Jacobi diffusion the Laplace transform 
of which can be computed explicitly in terms of Gauss' hypergeometric func- 
tions. It is also seen that similar techniques apply to compute the Laplace 
transform of an analogous functional for three-dimensional Bessel process 



The paper is organized so that in the next section some general results 
on the random time techniques connecting integral functionals and hitting 
times are presented. In this section we also discuss the Feynman-Kac ap- 
proach. Definitions and some properties of Gauss' hypergeometric functions 
are given in Section [HI In Section |31 we discuss some properties of Jacobi dif- 
fusions focusing on two special cases emerging from our computations. The 
functional in ()1.4p is treated in Section [HI and ()1.6|1 in Section El 

2 Transforming diffusions by changing time and 
space 

In this section we give a fairly straightforward generalization, needed in Sec- 
tionlH of the main result in Salminen and Yor [12] for quite general diffusions. 
We also discuss the corresponding Feynman-Kac approach. 

Let -B be a Brownian motion and consider a diffusion Y on an open 
interval / = (/, r) determined by the SDE 




(1.5) 



{Rt;t>0}: 




(1.6) 



dYt = cr{Yt) dBt + h{Yt) dt. 



(2.1) 
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It is assumed that a and h are continuous and a{x) > for x e /. The 
diffusion Y is considered up to its explosion (or life) time 

C := infO : Yt ^ I}. 

Proposition 2.1. Let g{x), x & I , be a twice continuously diiferentiable 
function such that g'{x) ^ 0. Consider the integral functional 

A:= f{g\YMYs)yds, te[0,C), 

and its inverse 

at := min [s : As> t], t e [0, A^). 
Then the process Z given by 

Zt:=g{Y,,), te[0,/lc), (2.2) 
is a diffusion satisfying the SDE 

dZt = dBt + G{g-\Zt)) dt, t e [0, A^). (2.3) 
where Bt is a Brownian motion and 

G{x) = {g'{x)a{x))-' Q aixf g'\x) + h{x) g\x)^ . (2.4) 

Proof. By Ito's formula for u < ( 

g{Y^)-g{Yo)^ g'{YMys)dB,+ {2-'g"{Y,)a''{Y,) + g'{Ys)b{Ys)) ds. 
Jo Jo 

Replacing u by at yields 

rat rat 

Zt-Zo= g'{Y,)a{Y,)dBs+ {g'{Ys)a{Ys)yG{Ys)ds. 
Jo Jo 

Since at is the inverse of ^4^ and A'^ = (^g'{Ys)a{Ys)Y we have 

< = ^ = i9'iyatHyat)r"- (2.5) 



From Levy's theorem it follows that 



Bt:= / g'{Y,)a{Y,)dBs, tG [0,71^). 



is a (stopped) Brownian motion. As a result we obtain for t < A, 



Zt-Zo = Bt+ [ {g'{Yja{YjYG{Yjda, 



= Bt+ [ G{g'\Z,))ds. 



□ 



It is interesting and useful to understand how the relation (j2.2|] expresses 
itself on the level of differential equations associated with the diffusions Y and 
Z. In particular, we consider in this framework the Feynman-Kac formula. 
First, we give the following proposition which is a generalization of the result 
for Brownian motion stated in [S| VI. 15 p. 113 and concerns distributions 
of functionals of Brownian motion stopped according to the inverse of an 
additive functional. 

Proposition 2.2. Let F and f be continuous functions on I. Assume also 
that F is bounded and f is non-negative. Let t be an exponentially (with 
parameter X) distributed random variable independent ofY. Then the func- 
tion 




is the unique bounded solution of the differential equation 



Y^u"{x) + b{x)u'{x) - {X{g'{x)a{x)Y + f{x))u{x) 

= -X{g'{x)a{x)YF{x). 



(2.6) 



Next, introduce for x G J := g{I) 




(2.7) 




^)Qix) 



(2.8) 
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By the classical Feynman-Kac result the equation (I2.8|l has only one bounded 
and non-negative solution, and, consequently, we have the following proba- 
bilistic interpretation 



Q{x) = ( F{g-\Z"^)) exp ( - f ^ 

where the diffusion Z° has the generator 

1 dx'^ dx 

and can be viewed as a solution of the SDE ()2.3p . Using ()2.7p and the 
probabilistic interpretations of U and Q we obtain 

^,^nV.)e.,{-[^J^,s)') (2.9, 

From this identity we can recover, e.g., that the distributions of and Z^ 
are equal (from ()2.2|) we, of course, know much more). In this paper, see 
Sections 4.2 and 4.3, the identity (I2.9jl is used in a particular case to find 
the solutions of the equation (j2.8j) in terms of the (known) solutions of the 
equation (j2.6jl . Notice also that Proposition 12.21 can be proved by using first 
the classical Feynman-Kac formula for ()2.8p and then Proposition 12.11 

Next two propositions connect occupation time functionals with first hit- 
ting times. We use the notation in Proposition 12.11 



Proposition 2.3. 1. Let x E I and y E I be such that P^-a.s. 

Hy{Y) := inf{t : = y} < oo. 
Then, under the assumptions in Proposition \2.1l 

Ah^^y) = inf{t : Zt = g{y)} =: Hg(^y){Z) a.s. (2.10) 
with Yq = X and Zq = g{x). 
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2. Assume g{r) := lim^^r g{z) exists (recall that r denotes the right hand 
side end point of I). Suppose also that each of the following statements holds 
a.s. 

(i) ( = +00, (a) lim Yt = r, (in) Aoo '■= lim At < 00. (2.11) 

t— >oo t— >+oo 

Then 

A^ = Hg^r){Z) a.s. (2.12) 

Proof. To prove (jSHni), notice that from the definition of Z (see (12. 2j) ) we 
obtain 

ZM=9{Yt). 

Letting here t — > Hy{Y) and using the fact that g is monotone yield ()2.in|l . 
The claim ()2.12j) is proved similarly. 

□ 

Remark 2.4. A sufficient condition for (iii) in (|2.1H1 is clearly that the mean 
of Aoo is finite: 

(Aoo) = _^ {{9\ys)o{Ys)f) ds 

Goix,y) {g'{y)a{y)fm{dx) < 00, 

where Go denotes the Green kernel of Y and m is the speed measure (for 
these see, e.g., Borodin and Salminen [3]). For a neccessary and sufficient 
condition in the case of a Brownian motion with drift, see (jLip . 

For the next proposition recall from p.5|l and p.2p that the inverse at of 
At is given for t < A^ hy 

at = f {g\g-\Z,)) a{g-\Z,)))-' ds. (2.13) 

Proposition 2.5. Assume that each of the following statements holds a.s. 
(z) C < 00, {ii) \miYt = I, (iii) = 00. 

t ^c^ 

Then 

{g'{g-\Z,))a{g-\Z,)))-' ds = C a.s. (2.14) 
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Proof. The claim is immediate from (|2.13jl using the assumptions (i)-(iii) 
and the fact that at is the inverse of At. Notice also that in this case ( = 
\im,^iH,{Y). □ 

We conclude this section by a proposition which relates the results in 
Propositions 12 . 31 and 12 . 51 to the Feynman-Kac formula. We consider only the 
case in Proposition 12.31 2 and leave the other cases to the reader. 

Proposition 2.6. 1. Under the assumptions of Proposition Wl\ 2. the func- 
tion 



^x) := E,. (^exp (^-p {g' {YMys)Y ds^"^ , 



X e I, p> 



is the unique increasing positive solution of the differential equation 

^ a'^{x) + h{x) ^'{x) - p {g'{x)a{x)Y <!}{x) = 0, (2.15) 

satisfying lim^..^^ ^{x) = 1. 
2. Define for x G g~^{I) 

^{x) ■.= ^g-\x)). 

Then ^ is the unique increasing positive solution of the differential equation 

^ ^"{x) + G{g-\x)) ^\x) ~ p^{x) = 

satisfying liiRx^r'^ix) = 1. Moreover, 

^{x)=E, {exp{-pHg^r){Z))) 

= E,-i(,) (^exp (^-p ^ {g'{YMys)f ds^ ^ . (2.16) 

Proof. The first part of the proposition can be proved as Proposition XX in 
Salminen and Yor |T2|. The fact that \l/ has the claimed properties is an easy 
exercise in differentiation. The first equality in ()2.16|) is standard diffusion 
theory and the second one follows from the definition of \E'. □ 
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3 Gauss' hypergeometric functions 

We start with by recalling from Abramowitz and Stegun [T| p. 556, the 
definition of Gauss' hypergeometric series (or functions): 



i^K/^,7,x). r(/3)r(a)^ r(7 + n) n\- ^'"'^ 



n=Q 



It is easily seen that F is well defined for |x| < 1 but we consider F only for 
real values on x such that < x < 1. However, we allow complex conjugate 
values for a, and /3, but take 7 real (see Theorems 15. II and 16. ip . Notice that 

F(a,/3,7;x) = F(/5,a,7;x). 

It is well known (see Abramowitz and Stegun [I] or Lebedev p. 162) that 
F is a solution of the ODE 

x{l - x) v"{x) + (7 - (a + /5 + 1) x) v'{x) — a (3 v{x) = 0. (3-2) 

Straightforward calculations show that also 

F{a, /3, 7; x) := F{a, /3, a + /3 + 1 - 7; 1 - x) (3.3) 

is a solution of (j3.2|) . 

Recall also that for < x < 1, and Re (7) > Re (a) > it holds 

F(a,A7;3:)= ^^/^ - ty-^-\l ~ tx)-^ dt. (3.4) 

r(a)r(7-a) Jo 

From ()3.4|) we obtain by changing variables 

F{a, l3, 7; x) = F{a, /5, a + /? + 1 — 7; 1 — x) 

= r^rr^t!"^! rt-^(l + t)^— ^(1 + tx)-/^c^t, (3.5) 

r(«)r(/5 + 1 - 7) Jo 

which is valid for < x < 1, Re(/5 + 1 — 7) > 0, and Re(a;) > 0. Notice that 

F(a,/3,7;1) = 1 

and 

^( R m r(a + /3 + i-7)r(i-7) 

^' ^' °) = r(a + i-7)r(/j + i-7) ' ^'-'^ 

where all the arguments in our case (see Theorem 15.111 are positive or have 
positive real parts. 
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4 Jacobi diffusion 



4.1 Definition and two particular cases 

Consider the diffusion X = {Xt : t >0} living in (0, 1) satisfying the SDE 

dXt = ^/2Xtil - Xt) dBt + i^-ia + l3+ l)Xt) dt. (4.1) 

Following Karlin and Taylor |H] p. 335, X is called a Jacobi diffusion with 
parameters a, /?, and 7. See the papers by Mazet ^j, Warren and Yor [T6] . 
Hu, Shi and Yor Schoutens [Hj, and Albanese and Kuznetsov for some 
results and applications of Jacobi diffusions (and also for further references) . 
Next proposition gives a basic simple property of Jacobi diffusions. 

Proposition 4.1. Let X be a Jacobi diffusion with parameters a, (3, and 7. 
Then Y := 1 — X is a Jacobi diffusion with parameters a, P, and a+/? + l— 7. 

The generator of X is given by 

Qu{x) = x(l — x) u"{x) + (7 — (« + /9 + 1) x)u\x), 

and we use 

S{x)= r y~\l-yr-^-^-'dy 

Jl/2 

as the scale, and 

m{xi, X2) = / y^'\l - yY+^-^ dy, < a;i < X2 < 1 

J Xl 

as the speed. Obviously, X exhibits very different behaviour when the values 
of the parameters a, P, and 7. are varied. In view of our applications we 
consider two cases. 

Case 1: 7 < 1 and a + P > 1. 

Here, 5'(0+) > —00 and 5'(1— ) = +00, and by the standard boundary 
point analysis (see, e.g., Borodin and Salminen j^) 

• if 7 < then is exit-entrance (regular) , 

• if < 7 < 1 then is exit-not-entrance. 
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• 1 is entrance-not-exit. 



Consequently, X is transient and 



Hq{X) := inf{t : Xt = 0} < oo a.s. 



In the regular case, i.e., 7 < 0, we choose to be a killing boundary. We 
remark that the boundaries and 1 are as displayed above also if the condi- 
tion a + P > 1 is extended to a + /3 > 7 — 1. 

Case 2: 7 = a + /? > 2. 

Now we have 5'(0+) = —00 and >S'(1— ) = +00, m(0, 1) < 00, and 

• is entrance-not-exit, 

• 1 is entrance-not-exit. 

Hence, X is positively recurrent. Notice, cf. Proposition 14. IL that if X^^^ 
satisfies (HH) with 7 = a + /? then X^^) := 1 - X^^^ solves (HH) with 7 = 1. 
This fact is also transparent in the expressions for S and m. 

4.2 Functionals of Jacobi diffusion; Case 1 

Consider the Jacobi diffusion given by ()4.H1 with 7 < 1 and a + P > 1. Let 



Clearly, / is decreasing, /(0+) = +00, /(I—) = —00, and 



We are interested in the process {/(Xj) : t > 0}. By Ito's formula for 



for c,e>0 




< X < 1. 



t < Ho{X) 




-V2 



1 - 7 + (g + /3 - 1)X, 
X,(l-X,) 



ds. (4.2) 



v/x,(i-x,) 
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Define for t < Hq{X) 

2 



At := — I ^, ds. 



and set At = +00 for t > Hq{X). Let a be the inverse of A, i.e., 

at := inf{s : Ag > t}, 
and notice that at < Hq{X) for all t > 0. Moreover, 



at = - I XaAl-XaJds. (4.3) 

The process given by 

9 cX„ 



Zt:=bog(^-^] (4.4) 



is well defined for all t such that at < Hq{X), and from (I4.2^1 we obtain 

^ ^ a + (3-l 



Zo = B^ + - 1-7+ "t^, . 1 ds, (4.5) 



where i?° is a Brownian motion. From ()4.5jl it is seen that Z is, in fact, 
non-exploading and, therefore, from (j4.4jl it follows that at < Hq{X) for all t 
and at Hq{X) as t ^ 00. Hence, it also holds that — > 00 as t — Hq{X), 
and we are in the case treated in Proposition 12.51 We have 

ceMOZ.) ,,^ri^(l^^u ds 



(c exp(^Z,) + 1)2 J, Xa^ V Xa, 



Xa(l-X,)ds 







2 



and the statement in Proposition 12.51 can be formulated as follows 

Proposition 4.2. Let Z, Zq = x, and X, Xq = l/(ce^^ + 1), be as above. 
Then 

9^ cexp(9Zs) ,,,,,,, , , 

JL (cexp(^l);i)^ '''-^°W (4.6) 
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As explained in SectionEl we can use the Feynman-Kac method to deduce 
that the Laplace-transforms of the functionals in ()4.6|1 are equal. However, 
because the present case is not covered by Proposition l2.6[ we formulate here 
a result connecting the solutions of the hypergeometric differential equation 
to the solutions of the equation induced by the generator of Z with the 
potential term as on the left hand side of (j4.6jl . 

Proposition 4.3. Let z{x) = 2;(q;, /5, 7, x), x G (0, 1), be an arbitrary solu- 
tion of the hypergeometric differential equation 

x{l - x)z"{x) + (7 - (a + /5 + l)x)z'{x) - apz{x) = 0. (4.7) 

Then for c > and 6' e R the function q{x) = z(y(x)) with 

y{x):={ce'^ + l)-' 

satisfies for x G (—00, 00) the equation 



9x 



Proof. It is sufficient to prove this statement for a = 1 since one can shift 
the argument. Differentiating the composition of the functions we have (the 
corresponding arguments are omitted) 

? = -v^ , / = 74^ + —r^ rr^^ • 



e 



+ 1)' (e^- + l)^ (e^- + l) 



Taking in (|4.7|) the argument y{x) instead of x, one obtains 
72^ + 7 j——]z' - a^z = ^. 



Finally, from these relations we have 

q [x] = — 2 ; 7 + -5 ]z + apz 



/e^" + a + /3 \ , ^^^/^e'^" , , 

= —6' s 7 9 a;) H o9 2;) 

V e'^^ + 1 ) ^ ' [^^Bx ^ x)2yv ; 

proving the equation (I4.8jl . □ 
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Remark 4.4. Notice also the converse of Proposition 14.31 if g is a solution 
of (j4.8jl then z{x) = q{y~^{x)), where 

2/-i(x):=^log(i— 0<a;<l, 

is a solution of (cf. ipT^ V 



4.3 Functionals of Jacobi diffusion; Case 2 

Assume that X is a Jacobi diffusion with 7 = a + /5 > 2 and recall that in 
this case X is recurrent. By Ito's formula 

log(l-X,)-l og(l--Y o) 



7-(. + /? + l)A-.^^,_l /'2.Y.(1-A-.)^^ 



1 - A'. 2 ;„ (1 - X, Y 



Let 



Q \i 1 — Xg Jq 1 — Xg Xg 



•= ^ / 



and define its inverse at := inf{s : Ag > t}. By continuity, At < 00 and 
at < 00 for all t and A* — 00 and a* — >■ 00 as t ^ 00; hence, we are in the 
case covered by Proposition 12. 31 1. Instead of simply refering to Proposition 
2.31 1. we give some details. Firstly 

«t = y / ^r^ds, (4.10) 



and the process given by 



Zt := -]-\og{l - XaJ (4.11) 
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is well defined for all t > 0. From (I4.9jl we obtain 

e f' a + (3 



Zt-Zo = B° + - / ds 

2 Jo exp{eZs)-l 

where B° is a Brownian motion. From (|4.1()|] and M.lljl we obtain 

t/ '^s 

For y > introduce 

HyiZ) := inf{t : Z^ = y}. 



Now, from (j4.1Hl . 

y = Zrr fy. = — - log (l — Xa 

and, hence, 

a^^(g)=inf{t : X, = l-e-'y}. 
Consequently, we arrive to the result (cf. Proposition 12. 31 1) 

Proposition 4.5. Let X be a Jacobi diffusion with •j = a + P > 2 and Z as 
defined in i)4.J jj) witli Zq = x. Tlien for y > Q 



o2 rHy{Z) ^ 1 

-j [expi9Z,)-lj ds = Hy,{X). (4.12) 



wliere y* = 1 — e ^ ^ and Xq 

It is useful (again) to give the corresponding result for differential equa- 
tions. Notice that in this result we do not interprete solutions probabilis- 
tically and, therefore, can formulate a more general statement without any 
restrictions on the values of the parameters. 

Proposition 4.6. Let z{x) = z{a, P,'j,x), x G (0, 1), be an arbitrary solu- 
tion of the hypergeometric differential equation \4.7\) . Then for 9 > the 
function s{x) = z [l — e~^^) satisfies for x e (0, oo) the equation 
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Proof. The proof is analogous to the proof of Proposition 14 .31 We have 
Taking in (j3.2jl the argument 1 — e instead of x yields 

- (« + /5 + 1 - Y^)'' - = 

Consequently, 

7 A , e^a(3e-^'' 



s"(x) = g'e-^"(a + /? + l- ~l)z' + - —z 



1 -e 



-Ox 



2^, n^-0x 



n( n 7 \ u \ a ape , , 

proving the claim. □ 
Remark 4.7. 1. We have also 
log TO- log (Xo) 



Choosing 7 = 1 and letting 



6' k X, 



it is seen that the process 



where a° is the inverse of A°, satisfies the same SDE as Z. 
2. Let z{x) = 7,x), x G (0,1), be an arbitrary solution of the 

hypergeometric differential equation ()4.7|) . Then for 6 > the function 
q{x) = ^(e~^^) satisfies for x > the equation 

^ (^) + ^(l - 7 + ^ )g^(^) - = 0. (4.13) 
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5 Perpetual integral functional of BM(//) 

Our first main result completes, in a sense, the result in jTTj, see also [12], con- 
cerning the translated Dufresne's functional but, moreover, it gives Laplace 
transforms for many new perpetual integral functionals. The functional we 
analyze is 

/oo(p,g):=^ ( ^ exp(^^i?^)) + 1 ^ (c exp(^^5i^)) + l)') 

The notation It{p, q) is used when the integration is from to t and lHy{p, q) 
when t equals 

Hy := inf{s : = y}, 

the first hitting time of y. 

Theorem 5.1. Let B^^^^ be a Brownian motion with drift /i > started from 
X. Then for c>0, 9>0, p>0 and p + q >0 

E^(^exp - /oo(p,g))) 

= K v{x)'' F{a, (3,a + f3 + 2^i/9 ; v{x)) (5.1) 
where F is Gauss' hypergeometric function as defined in and 

c exp{9 x) 



k= {a + (3~l)/2, v{x) 



c exp{6 x) + 1 



T{a + 2f,/9)T{P + 2fi/9) 

T{a + (3 + 2fi/9)T{2fi/9) ^ ' 



and 



a = --fi/9+y/fi^ + 2{p + q)/9 + 2 V 1 + 8?/^'> (5-3) 

1 ,„ ; r,„ 1 



P=--^/9+y/^,^ + 2{p + q)/9 - 2 V 1 + 8g/^'- (5-4) 



Proof. Recall, e.g., from [12] that 
^{x) := E^(^exp - I^{p,q) 



E,(exp (^p[l- v{b[^^)) +q(l- v{b[^^) \ \ dt 
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is the unique positive bounded function such that 

^^"(x) - (^p{l- v{x)) +q{l- v{x)f^ ^(a;) = (5.5) 

and hm^-^+oo "^{x) = 1. We remark that in [T^ it is required \l/ to be increas- 
ing, but from the proof therein it is clear that an equivalent requirement 
is boundedness. To find \1/ we use Girsanov's theorem and Proposition 14.21 
Firstly, recall that the process Z as defined in (14. 4|] solves the SDE 

dZt = + ^(l - 7 + (« + /? - 1)(1 - v{Zt))') dt. 

Choose now 7 such that /i = 6{l—'-^)/2. By Girsanov's theorem, the measures 
and P^'"-* induced by Z and B^^\ respectively, and defined in the space 
of continuous functions uj : R+ ^ R are locally absolutely continuous with 
the Radon-Nikodym derivative given by 

P^l^, = exp(A)P(^M^., (5.6) 

where Tt is the cx-algebra generated by the coordinate process up to time t 
and 



with 
Observe 



A:=^ /K)rf^,-i^ (J\uJs) + 2^lf{uJs))ds (5.7) 



r 1 r^'iy) 1 

V{x):= j {I - v{y)) dy = - j — ^ dy =- log v{x). 

For the first integral in ()5.7|1 we obtain by Ito's formula 

1 /■* 

(1 - v{ujs)) dus = V{uJt) - V{uJo) + - v'{ujs) ds 

2 Jo 

9 /"* 

= V{uJt)-V{uJo) + - v{ujs)il - vii^s))ds 



2 _ 

V{uJt)-V{uJo)-^ 1 {l-v{uj,)fds 



2 
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+ - / {1 -v{ujs))ds. 



Letting (cf. Proposition I4.2jl 



2 



(c exp(6'co's) + 1) 



ds 



we have 



It = —I v{ujs){l - v{ujs)) ds 



= y - viujs)) ds-^^ v{ujs)y ds. 

The absolute continuity relation ()5.6p yields 



E 



f(exp(-r/,)) =E('^)(exp(A-r/i 



(5.8) 



Straightforward but lengthy computations show that choosing r = a/5 the 
claimed expressions llb.'M and (I5.4jl for a and /3, respectively, are such that 



exp ( ek V{Bi^^) - ek V{x) ) exp ( - It{p, q) ) . (5.9) 



E, 



We remark that r = a/? can attain both positive and negative values. The 
equalities (|5.8j] and ()5.9|1 hold also for the first hitting time Hy with y > x, 
i.e., 



Ef exp(-r/H, 



vix] 



exp - iHyip^q) 



Letting y +oo (remember that Zt 



-oo and B 



(m) 



we obtain by monotone convergence and Proposition 14.21 

$(x) = {v{x))^E^J exp ( - a (3 Ho{X) 



+00 as t — > +cxd) 



(5.10) 



where x* := 1 — f (x). Recall that ^/'(x) ^ 1 as x ^ +oo. This can also be 
verified from (j5.1()j) by observing that as x — > +oo then x* — >^ and the right 
hand side of (j5.1()jl tends to 1 because is an exit (or killing) boundary point 
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for X. Notice also that if p + q = (see Example 5.3) (in this case a(3 > 0) 
then k = and we have 



< lim "^ix) 



lim Ef. 



exp ( -apHo{X) 



Also if p + g > then \E'(x) — > as x ^ — oo. In this case, because v{x) 
we cannot without further analysis claim that 



lim E^. 



exp ( -aPHoiX) 



< +00 



in case a/3 < (but this will follow from our discussion !). For the general 
theory of diffusions we know that the function 



^^^(x*) ■.= E^.[expi^-a(3Ho{X) 

is a solution of the hypergeometric differential equation 

(Clearly, ip^^ as a function of x* is decreasing if aP > 0, and increasing 
if aP < 0. Hence, because 1 — v{x) is decreasing, ip^^ as a function of x is 
increasing if a/3 > 0, and decreasing if a/? < 0. In the latter case, multiplying 
ipa/3 with the increasing function v{x)^ makes the product increasing.) 

Let F be the hypergeometric function as defined in ()3.H] with a and P 
as in (I5.3^1 and ()o.4jl . respectively, and 7 = 1 — {2fi/9). Notice that a and P 
can be conjugate complex numbers and 

a + P > 1, and 7 < 1. 

Moreover, as is easily seen, 

Re{p + 1 - 7) = Re{p + {2^/9)) > 0. 

From Section El it now follows (cf. (|3.5|) ) that the function 

X F(a, P,a + P + l--f;l-x) = F{a, P,a + p+ (2/i/6'); 1 - x) 

is a bounded non-negative solution of the hypergeometric differential equa- 
tion. Proposition 14.31 and formula yield 



^°(xl 



{v{x))^ F{a, p,a + p + {2fi/9); v{x)) 



is a nonnegative bounded solution of (|5.5p . Consequently, by the uniqueness 
we have 

^(x) = ^°(x)/^°(+oo), 
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where, from 

^°(+oo) := lim 



X—r + OO 



= F{a,(3,a + P+i2fi/e);l) 

_ r{a + (3 + {2fi/e))r{2fi/e) 

~ r(a + (2/i/^))r(/3+(2/i/e))' 
as claimed. □ 
Example 5.2. Choosing in 1^5.1^ p = 0, q = ■y /a^ , 9 = 1, and c = 1/a gives 

E,(exp(-7^ (a + exp(5(^)))-2rfs)) 

= K v{x) (V«V2+27-a/^)/a ^(^t^ « + ^ + 2/i ; t;(a;)), 

where 



1 / „ 27 /I 27 ^ 1 / „ 27 /I 27 

"=2-^+V^ +^ + V4 + ^' ^=2-^+V^ +^"V4 + ^' 

-ftT is given by (j5.2jl with a and /3 as above, and 

exp(x) 



v{x) 



a + exp(xj 



Notice that if /x = 1/2 then /5 = and because F{a, 0, 7; x) = 1 for all \x\ < 1 
(see (j;i2|) ) we obtain the result in [T2j (see also [TT]): 

E,(exp(-7 J {a + exp{Bi'/^y)r'ds)) = ^(3,)(2a)- (V^-). 

From Proposition 12.31 2 it follows that 

POO 

/ {a + exp{Bl'''^))-^ds = HoiZ) a.s. 
where Z satisfies the SDE 

/ 1 ae"-^* \ 
dZt = dBt+ ( (/i - -) ^ _ ^^zt + ] 

The function g is in this case 

g(x) := -- In (1 + ae"^). 
a 
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Example 5.3. Let in (jSIIl) p + g = 0, p = 4p > 0, 6* = 2, and c = 1. Then 



1 1 . . 1 1 



r(/x)r(/i + i) ' 

and we obtain the result (see Vagurina [T5] ) 



By Proposition 12312 

/•°° ds 
Jo cosh\Bi''^) 

with Z determined via the SDE 



H^iZ), a.s, (5.11) 



dZt = dBt + ( - ctn Zt H — . ^ ] (it, = 2arctan exp(i?Q'^^). 



^ 2 sin 
In this example 



(/(x) := 2arctane^', 5''(a;) 



cosh X ' 
and, therefore, 

^/ -1/ 1 A* 1 1\ a; 1 / 1\ X 

Giq (X)) = - ctnx H = - u tan — I — u H — ctn — . 

^ " 2 sinx 2 y 2j 2 2\ 2) 2 

Further, 

E,(exp(-pif,(Z))) = KF{a, /3, 1 + p; sin2(^/2)) 

Example 5.4. Take p = he and g = a^c^/2 in the defintion of Ioo{p,q). 
Then, letting c — > +oo, we obtain by monotone convergence (using JUII)) 

/■°° / he (? \ 

™ cexp(^B^^) + l ^ ^ (cexp(^5i'^)) + l)'^ 



/ /i exp(-e5('^)) + — exp {-2eB^f^)\j ds 
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The Laplace transform of the functional on the right hand side is given in [2] 
formula 2.1.30.3(2) p. 292 and can written as 



exp [h exp(-05(^)) + y exp(-2e5(^))^ ds 

_r(l/2 + ,/. + Va^)^^ ,_c 



E, 



r(2/i/0) V e 

where the Kummer function U is connected to the Whittaker function W via 

W^^mix) = W„,-^(x) = g-./2 _ ^ ^ ^/2, 1 - 2m, x), 

see Abramowitz and Stegun [l]. To obtain the formula ()5.12|l from (|5.H1 
observe first that as c ^ +00 we have a ~ 2ac/9, 

^ 2 ^ a^' 

and 



'-^^^J -exp(^--e- 
and recall that for Re7 < 1 and Re{a + P — •y) > —1 

Fr^^T^TT^^("' /5, « + /3 + 1 - 7; 1 - -) = f/(/?, 7, x), 
a^oo r(a + p + 1 — 7) 

see Abramowitz and Stegun p. 

6 Perpetual integral functional of BES(3) 

Let R = {Rt : t > 0} denote a 3-dimensional Bessel process (or, equivalently, 
of index 1/2). The generator of R is 

1 d'^u ^ 1 du 

2 dx'^ X dx 

In the next theorem we give the Laplace transform of the functional. 

(exp(^^.)-l + (exp(J,)-l)0^'- 
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Theorem 6.1. Let R be a 3-dimensional Bessel process started from x > 0. 
Then for 6 > 0, p > and q > 



F{a, (3,^,1 -e-'^) (6.1) 



Ex exp - /oo(p,g) 

(l_e-^-)7/2r(a)r(/3) 



(l_e-^-)7/2 /•! 



t''-\l-tf'\l-tx)''-^dt, 



where 



V^Vl + 8g/^^2+ 1^2(g-p) (6.2) 



a 



2 2" 9 



f3 = ^ + ^^l + 8q/e^-^,/2{^), (6.3) 



and 7 = a + /5 = 1 + VI + Sq/O^. 

Proof. To start with, let Z° be a diffusion witfi tfie generator 
wliere 7 > a + /? > and 

Notice that if 7 = a + /? then Z° coincides with the diffusion Z introduced 
in Section 1131 The measures induced by Z° and by R in the canonical space 
of continuous functions are absolutely continuous with respect to each other 
when restricted to the cx-algebra J-'t generated by the co-ordinate mappings 
up to a fixed but arbitrary time t. Let Pf and be the measures associated 
with Z° and R, respectively, when both processes are started from x > 0. 
Then 

Pf 1^, =exp(A)Pf (6.4) 
where the exponent Dt in the Radon-Nikodym derivative is given by 

Dt = - [ (uj;^ - h{us)) dujs + l [ (cu;2 - h\uj,)) ds (6.5) 
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Consider the stochastic integral term in Under the measure the 

co-ordinate process is the Bessel process R started from x. Recalling that the 
quadratic variation of R is t, we obtain by Ito's formula 

(i) - £ R-' dR, = log (|-) - i £ Rf ds 
(") £ KR,) dR, = 1 log (^^) - ^-^^^{Rt - X) 

Next consider the bounded variation part in We have 
(m) -- / h\Rs)ds = / —J- —rds ' ' t 



BRs 

ds. 



2Jo ' ' 8 7o (e^^— 1) 



Because Rt — > +oo as t — > +oo, it is seen from (ii) and (iii) above that taking 
■J = a + j3 will lead to a remarkably simple special case. Indeed, if 7 = a + /3 



In view of Proposition 14.51 we consider now the functional (assume that 
the continuous function u is such that the functional is well-defined) 



% ■=— {exp{eujs) - ds. 
^ Jo 

By absolute continuity, 

Ef (expi-rZ)) = (exp(A - rl, 



t) , 
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and, further, for Hy := mi{t : Uf = y} with y > x 

Ef (exp(-r/Hj) = Ef {expiDn, - rlu,] 

Proposition 14.51 gives now 

Ef. {eM-rHy,)) = Ef (expp^, - r/^j) , (6.6) 

where x* = 1 — e~^^, y* = 1 — e^^^, and X is a Jacobi diffusion with 
parameters a, P, and 7 = a + /9 > 2. The identity ()6.6p is equivalent with 

Ef (exp(-/H,(p,g))) = ^ Q~^„'y ^ Ef (exp(-rif,0) , (6-7) 
where 

02 q2 

P=j(2?^-7) and g=— 7(7-2). (6.8) 

Letting r = ajd it is seen from (|6.8jl . after straightforward computations, that 
a and /5 can be expressed as in (16. 2j) and (j6.3|) . respectively. To conclude the 
proof it remains to compute the Laplace transform on the right hand side of 
()6.7j) and let y +00. Because x < y we have also x* < y*, and. hence, 

Ef (exp(-a/5 //,.)) = m^y, (6-9) 

where ipai3 the unique (up to a multiple), positive increasing solution of 
(13. 2j) with J = a + p. We remark that the uniqueness follows from the fact 
that both boundaries, in this case, are entrance-not-exit. Moreover, because 
Hy* — > +00 as ?/* — s> 1 (also as y* 0) it must hold that 'ipa/siv*) ~^ +^ 
as 7/* ^ 1. Similarly, there exists an unique positive decreasing solution ip^j^ 
such that ^Pa/siy*) ~^ +^ as y* 0. From the general theory of differential 
equations it follows that all other solutions of (i;i2|) can be expressed as 
linear combinations of ^/^^ and y^f^. Recall (see Section ISJ that the function 
X I— > F{a, x) is a solution of ()3.2p . Now 7 = a + and from Abramowitz 
and Stegun [T] 15.3.10 p. 559 

F(a,/3,a + /3;x) ~ -^y^^-^-^ log(l -x), as x T 1, (6.10) 

1 (a)l [p) 
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and, because F{a, P,a + P;0) = 1, it follows that x i— F{a, P,a + P; x) 
increasing, and, consequently. 

Therefore, from ()6.9jl and ()f).7|] . 

T _y fl-e-'^y^' F{a,P,a + P;x*) 



Letting here y — > +00 and using (|6.inil proves the claim. 
Remark 6.2. The decreasing solution ip-^^^ is given by 

^api^) = -^(«) (3,a + /3; x) := /?, 1; 1 - x). 

From the definition of F and the facts that, in this case, a/? > and 

(a + k){(3 + k) = a(3 + k{a + (3) + A;^ 

it is seen that F has the desired properties, i.e., 

lim /3, a + /5; x) = +00 and /3, a + /3; 1) = 1, 

x— >o 

implying that F is decreasing. 

Example 6.3. In the particular case p = q = 4:p, and 6* = 2, we have 

1 1 



« = /?=- + - Vl + 8p 



and, therefore 



Example 6.4. Choosing g = yields 



a = 1 + ^ V^, /? = 1 - i =: «. 
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Hence, a + (3 = 2 and 



1 _ „-ex 

— T(a) T(a) F(a,a,2;l- e"^^) 

1 - e"^^' n^/2p 



where the formula (see P formulae 6.1.28 and 6.1.29. 

r(i + iy)r(i-iy) 



F(a,a,2;l-e-''^), 



sinh('7r?/) 
is used. Letting x — > yields 

Eoexpf- f ds) = . . w"^... - (6.11) 



exp{eRs) - 1 / e smh{n^/2^/9) ' 

The term on the right hand side is the Laplace transform of H^^/g^R) where 
the BES(3) process R is started at 0, i.e., ()6.1ip is equivalent with the fol- 
lowing identity due to Donati-Martin and Yor ||4j p. 1044: 

Jo exp{eRs) - 1 

The derivation of (I6.12jl in jlj is very different than the one presented above, 
and is based on a formula for the Laplace transform of an integral functional 
of a two-dimensional Bessel process. 
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